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We offer a method of correlations mapping on the full celestial sphere that allows to check the 
quality of reconstructed maps, their non-Gaussianity and conduct experiments in various fre- 
quency ranges. The method was evaluated on the WMAP data, both on the reconstructed maps 
and foreground components, and on the NRAO VLA Sky Survey (NVSS) data. We detected a 
significant shift in the correlation data of the dust component, which can be preconditioned by a 
more complex dust model than the one currently in use for component separation. While study- 
ing the NVSS correlation data, we demonstrated that the statistics of the coinciding spots in the 
microwave background and in the NVSS survey corresponds to the one expected in the ACDM 
model. This can testify for a chance coincidence of the spots in the NVSS and WMAP data in 
the CMB Cold Spot region. Our method is software-implemented in the GLESP package. 
PACS: 95.75.-z, 98.62.Ve, 98.70.Dk, 98.70.Vc, 98.80.-k 
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1. INTRODUCTION 

Full celestial sphere surveys in different frequency 
ranges have led to a major advance in the research 
of global properties of the world around us. Such im- 
portant matters as determining the cosmological pa- 
rameters of the Universe in its early stages, the for- 
mation of the large-scale structure of the Universe 
and the phenomenon of dark energy are among them. 
The observations of the cosmic microwave back- 
ground radiation (CMB) by the Wilkinson Microwave 
Anisotropy Probed (WMAP) [D El [3], were revolu- 
tionary in modern cosmology. The data were recorded 
in five bands: 23 GHz (the K-band), 33 GHz (the Ka- 
band), 41 GHz (the Q-band), 61 GHz (the V-band) 
and 94 GHz (the-W band) with the measurements of 
intensity and polarization. The result of the analy- 
sis of these data were the CMB maps of anisotropy 
and polarization, the maps of foreground components 
(synchrotron and free- free emission, dust radiation), 
their power spectrum were as well deduced. For the 
harmonics of not very high order {£ < 100), a map of 
CMB anisotropy distribution is listed, reconstructed 
from the multifrequency observations of foreground 
components implementing the Internal Linear Com- 
bination (ILC) method [Tj. 

In present-day full sphere observations, the in- 
frared [4], radio [5] and optical [6^ wavelength ranges 



^ |http : //lambda.gsf c .nasa. gov] 



survey maps are actively used along with the mi- 
crowave background data. With the aid of correlation 
between the CMB and radio emission, one can study 
both the contribution of radio sources into the com- 
mon background [71 [5] , and the correlation properties 
of the radio source and CMB distribution occurring 
both on the large scale {9 > 2°), like the Sachs- Wolfe 
effect [9], and on the small scale (mainly 6 < 4') [TUj . 
like the Sunyaev-Zel'dovich effect. Correlations of in- 
frared and optical surveys (specifically with known 
red shifts) with CMB data allow us to study the for- 
mation of a large-scale structure. 

Another important application of correlation 
methods are studies of CMB maps quality and their 
cleaning level from other kinds of radiation. For ex- 
ample, the presence of correlation between the CMB 
and foreground components I12j carried in by the 
Galaxy, points to a residual contribution of the inter- 
fering emission in the procedure of components sepa- 
ration. Therewith, the presence of correlated compo- 
nents in a number of multipole frequency ranges leads 
to a change in the CMB signal statistics, that mani- 
fests itself through non-Gaussianity [1^. That, in its 
turn, makes the analysis of the maps of the studied 
signal and their power spectrum more complex. Quite 
a significant number of papers has been dedicated to 
the studies of statistic properties of the signal on the 
ILC map and to the discussion of its non-Gaussianity. 
These works were conducted using different methods, 
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like the phase analysis method [121 HSl E], Maxwell's 
niultipole vectors method [15], wavelet analysis [El 
[TTl HI] and Minkowski functional [19l|20]. In the 
non-Gaussianity studies, the correlation method al- 
lows separating a non-Gaussianity of a certain type, 
e.g. the one preconditioned by systematization during 
data processing. 

In this paper we generalize the approach of map 
correlation in various angular scales, earlier described 
in [211 122] ■ We offer a method of tessellated mapping 
and maps correlation visualisation on the full sphere 
for different angular scales of correlation search. We 
use this approach for the construction and study of 
the ILC and WMAP5 foreground components corre- 
lation maps, where the WMAP data are from the fifth 
year of observations [2]; as well as for the correlation 
maps of ILC and NRAO VLA Sky Survey (NVSS) 
radio sources [5]. 

2. THE METHOD OF CORRELA- 
TIONS MAPPING 

The correlation of two maps on the sphere is described 
by a correlation coefficient for a multipole i as follows 

I 

m = i,^^ J , (1) 

m—~£ m——£ 

where tem and Sim are the variations of CMB tem- 
perature and some other signal in a harmonic rep- 
resentation, * is a symbol of conjugation, £ and m 
are the numbers of spherical harmonics (of the multi- 
pole) and its modes in a harmonic decomposition of 
the signal on the sphere: 

X] H +a^,-my«,-m(6',0)) , (2) 

t m=l 

where Y^,m are the spherical functions, (0, 4>) are 
the spherical coordinates, ai^m are the coefficients of 
spherical harmonics satisfying the following condition 

y,,_™(0,</))-(-irr,:„(0,0), 

a,,™ = (~ir4„,„. (3) 

The value of the K{£) coefficient allows to check the 
harmonics correlation on the sphere, i.e. to compare 
the properties of maps in a given angular scale. How- 
ever, while looking for the correlated regions which do 
not recur in other regions of the sphere, this approach 
brushes out such single areas while common averag- 
ing over the sphere within one harmonic. In this case 



it is practically impossible to separate the correlated 
regions. 

We offer an approach realized within the frame- 
work of the GLESP package [13 [24], which allows 
detecting the correlations in the studied areas of a 
given angular scale for the maps, pixelized with higher 
resolution. The method is implemented in the pixel 
parametric space. The correlation mapping consists 
of assigning the pixel number p with a correlation 
result inside the solid angle Sp, calculated for two 
maps having higher resolution. As a result, we obtain 
a new map, in which the value assigned to each pixel 
reflects the correlation level of the studied maps in a 
given area. 

The correlation coefficient of anisotropy of the 
CMB temperature and a given signal for every pixel 
p [p — l,2,...,A'o, where Nq is the total number of 
pixels on the sphere), corresponding to the solid an- 
gle 5p and deduced for the maps on the sphere with 
an initial resolution, determined £max (hereafter the 
value of £max = 150 is taken for the original maps), 
there are 

E E (Ar(e„ M - Ar(Hp))(s(e., - Wp)) 



where AT{9i, (f>j) is the anisotropy value of the CMB 
temperature in a pixel with the coordinates {Oi,(pj) 
for the given pixelization resolution of the sphere, 
S{0i,(j)j) is the value of another signal in the same 
area, Ar(Sp) and S'(Sp) are common values in the 
area Sp, deduced on the data from the maps with a 
higher resolution, plotted by £max, o'ATp and as^, are 
the corresponding standards for this area. 

3. CORRELATION MAPS OF ILC 
WMAP5 AND FOREGROUND 
COMPONENTS 

One of the applications of the given approach is check- 
ing the cleaning quality of the CMB signal isolated 
from multifrequency observations. The presence of a 
residual (correlated) signal in the CMB data may lead 
to a bias in power estimation in a number of multi- 
pole ranges [5S| , which in its turn leads to a decrease 
in accuracy while determining the cosmological pa- 
rameters. 

Using the WMAP data — a map of synchrotron 
radiation in the K-band, a map of free-free emission 
in the V-band, dust radiation in the W-band and, 
finally, the ILC map — we have built correlation maps 
of foreground radiations and ILC at different angular 
scales. 
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Figure 1: Correlation maps of ILC and galactic foregrounds: dust, free-free and synchrotron radiation (from 
left to right), for different correlation windows: 16^ , 324' ^.''^d 5^0 (from top to bottom). The tone of the pixels 
reflects the value of the correlation coefficients from -1.0 to 1.0. The brighter the tone the smaller the absolute 
value of the correlation coefficient. 



As it was demonstrated earlier in the literature 
dedicated to the analysis of the ILC map and its sig- 
nal statistics [TTl [121 HSj j there are some serious ar- 
guments for the fact that in this map at different an- 
gular scales there exists some residual contribution of 
foreground components, which, among other factors, 
gives the found non-Gaussianity. This contribution 
may as well manifest itself in the earlier discovered 
connection in the quadrupole between a cleaned map 
of the microwave background and galactic radiation 
components [25l |26l [27| ■ According to the papers of 
the WMAP team [TJ H] , the ILC map is not intended 
for the CMB study on high multipoles, but it could 
be used for a foreground components analysis. This is 
the map we use for the foreground properties analysis. 

In our previous work 21j, using one-dimensional 
sections, we employed a WMAP mission foreground 
maps correlation search approach, and discovered in 
particular, that in the WMAP maps sections on dec- 
lination 5 = 41° there is a signal, correlated and 
anti-correlated with the data of the components be- 
ing separated. Among other things, on the angular 
scales indicative for the Galactic plane (multipole 
range £=10-20) some notable, compared to the simu- 
lated maps, correlations with the galactic foreground 
components were discovered. These correlations may 
be the consequence of the residual errors at the signal 
separation phase. 



In the present work we constructed the ILC and 
WMAP5 foreground components correlation maps on 
the full sphere. We as well estimated here the sta- 
tistical value of the isolated correlated areas using 
the Monte-Carlo method, comparing them with the 
Gaussian random fields maps, generated on full celes- 
tial sphere for the ACDM cosmological model [3] . The 
results of correlation mapping of the ILC and galac- 
tic foreground maps are demonstrated in Fig. 1 for the 
correlation windows with the sides of 162', 324', 540' 
(consequently, imaxi — 33, 16, and 10) for the three 
phase components. The window dimensions were cho- 
sen in such a manner that they fit into the range of 
angular scales of the Galaxy's influence (£=10-20), 
and the possible Sachs- Wolfe effect {£ < 45). 

The Galactic plane can clearly be seen on the 
maps presented in Fig.l. An important marker here 
is the correlation value statistics for every pixel. We 
constructed the histograms for all the given maps 
and compared them with analogous results for the 
simulated maps. The simulated maps were generated 
with an assumption that the CMB data are the re- 
sult of a Gaussian random process with an angular 
power spectrum, corresponding to the cosmological 
ACDM model. In total, there were constructed 50 
random simulated maps using the cl2map procedure 
of the GLESP package [13 [IJ, and with their aid 
we were able to evaluate the admissible interval, de- 
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termined by the correlation values spread dispersion. 
The pixel value distribution is presented in Fig. 2. The 
admissible limits of the correlation level constructed 
on the modelling data, are shown with dashed lines. 

Two foreground components correlations, dust 
and synchrotron, exhibit a shift in the range of anti- 
correlations. While in case of the dust component 
this shift is observed in all the studied scales, in 
case of synchrotron radiation it is presented on the 
histograms for 324' and 540'. The free- free emission 
histogram shows the distribution within the normal 
range, corresponding to CMB as a Gaussian ran- 
dom process. We have to note, that an analogous 
shift in the distribution was detected for low multi- 
poles at the analysis of the ILC signal reconstruction 
method [551 IHl via regenerating the simulated maps 
with a further intermixing with the foreground com- 
ponents, reconstructing the signal, then correlating 
the reconstructed maps and constructing the corre- 
lation coefficients distribution. In this approach the 
distribution evaluation consisted of two steps only: 
the correlation map construction and the calculation 
of the pixel values distribution histogram. 

A shift in the correlation coefficient distribution 
may be linked to an underestimation in the model 
used of the dust distribution in our Galaxy, where the 
excess in the negative coefficient values supports the 
hypothesis that in the corresponding pixels the cor- 
related CMB and dust signals have an opposite sign. 
This fact can be used for the correction of dust com- 
ponent contribution. We are planning further studies 
of the detected bias. 

4. ILC WMAP5 AND NVSS CORRE- 
LATION MAPS 

Let us consider another correlation mapping applica- 
tion. Among the CMB map research works there are 
those dedicated to the dark matter studies using the 
Sachs- Wolfe integral effect [9] , registered while corre- 
lating the CMB maps and the NVSS source distribu- 
tion (28^29]. 

The NVSS survey [5] is the most comprehen- 
sive Northern sky survey. It was conducted using the 
NRAO Very Large Array telescope (VLA) at a fre- 
quency of 1.4 GHz from 1993 to 1996, it has high 
sensitivity and covers the sky north of the declination 
5 = -40° (338840° or 82% of the celestial sphere). 
This survey is actively used for various statistical 
studies in cosmology. The catalogue of this survey 
contains 1.8 x 10^ sources and, according to it descrip- 
tion, it is 99% complete up to the integrated flux den- 
sities over S'i.4ghz = 3.5 mJy and 50% complete up to 
the flux densities of 2.5 mJy. The survey used the D- 



configuration of the VLA radiotelescope, and the size 
of the synthesized direction diagram in the half-power 
level, determining the resolution, constituted approx- 
imately 45". The survey data are available from the 
NRAO websitel, the virtual telescope SkyVie-wH and 
from the CATS databas«Q [30l [31] . 

Some non-Gaussian peculiarities were detected 
from the WMAP and NVSS data, poorly matching 
with the inflation ACDM model. Among them, for 
example, is the Cold Spot (CS), having the size of 
about 10° and the galactic coordinates of (/=209°, 
6=-57°). The CS, being statistically isolated in the 
CMB distribution, is one of the maps' particularities 
that contradict the hypothesis of homogeneous Gaus- 
sian background fluctuations. Originally, it was noted 
as a deviation from the Gaussian statistics while us- 
ing the wavelet analysis 46., _18j for the flrst year data 
recorded by the WMAP mission. Subsequently, in the 
NVSS maps there was detected a significant dip of 
the radio source density of matter in the CS region 
131], that allowed to hypothesize the existence of a 
Great Void 140 Mpc across, at redshift of z < 1, 
causing a gravitational anomaly leading to the inte- 
grated Sachs- Wolfe effect [9], and manifesting itself 
as a CS. At the same time there appeared some se- 
rious facts arguing that in this map at different an- 
gular scales there is some residual foreground compo- 
nents contribution, which produces the detected non- 
Gaussianity. This contribution may manifest itself in 
the earlier discovered link in the quadrupole between 
the cleaned CMB map and the galactic radiation com- 
ponents [221 [57] and, for another hand, determine the 
low multipole properties £ < 20, resulting in their un- 
stable reconstruction [25l[27j. Specifically, the CS par- 
ticularities may be explained via the effects of these 
multipoles: the deviations of peak clusters statistics 
around the spot [33], namely, an increase in the the 
positive peaks number. An independent study of the 
properties of a spot, detected in the region with close 
coordinates on the maps and the radio source count 
in the NVSS survey [5], has demonstrated that the 
studied Cold Spot, with its gigantic size and its mere 
existence being difficult to explain in terms of the ex- 
isting cosmological ACDM model, may as well be a 
simple statistical deviation, caused by the systematic 
effects |34j. 

Earlier we made a correlation analysis of the 
NVSS and WMAP 5 surveys, preliminarily modify- 
ing the NVSS map making use of the method of 
one-dimensional sections for different angular scales 



^ |http : //www . c V ■ nr ao . edu/nvss/ | 
^ littp : //skyview. g sfc ■ nasa.gov] 
^ |http : //cats ■ sao ■ ru] 
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Figure 2: Distribution of correlation index values in the pixels of the ILC and galactic foregrounds correlation 
maps: the dust radiation, the free-free and synchrotron emission (from left to right) for different correlation 
windows: 162! , 324' '^^'^ 54^0 (from top to bottom). The admissible limits of correlation level, constructed on 
the modelling data, are shown with dashed lines. 




Figure 3: ILC and NVSS correlation maps for different correlation windows: 162! , 324' "^^c^ 540! (from right 
to left). In the right bottom part of each figure there is an screened area, which was not considered in our 
analysis due to the lack of NVSS observations in the 6 < —40° range. 



[22j . To construct the modified NVSS map, we used 
a characteristic, represented by a mean square of the 
source flux density in an area of given dimensions 
with the center at the pixel, determined by the se- 
lected area. As a result we discovered that the corre- 
lation statistics in the studied scales (0.75°, 3°, 4.5°, 



and 6.75°) does not vary from the expected random 
maps. Therefore, we may conclude that the found cor- 
relations may be caused by a statistical coincidence. 

In the present work we carried out the map- 
ping of the correlation properties on the sphere us- 
ing the same data from WMAP5 and NVSS. The re- 
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Figure 4: Distribution of correlation index values in the pixels of the ILC and NVSS correlation maps for 
different correlation windows: 162! , 324' j '^i^d SJ^ff (from right to left). The admissible limits of the correlation 
level, constructed using the modelling data, are shown by the dashed lines. 



suits are shown on Fig. 3 and 4. While analyzing the 
map statistics (Fig. 4) we were taking into account 
the absence of NVSS data in the declination area of 
(5 < —40°. The NVSS map was constructed using two 
methods: (1) via computing the number of sources in 
every pixel and (2) via summing up the flux densi- 
ties of the radio sources, contained in the pixel area. 
As the correlation properties are practically undis- 
tinguishable (the difference of less than 1%), we shall 
demonstrate the whole analysis for the first case. 

As we can see from Fig. 4, the statistics of the cor- 
related pixel values lies within the admissible random 
deviations, what conflrms the earlier estimations |22j . 
This implies that in the correlation coefficient distri- 
bution no signal is found with a level higher than la 
on the given scales, caused by, for example, an in- 
teraction with the large-scale structure and dark en- 
ergy, manifesting themselves through the Sachs- Wolfe 
effect. On the other hand, as the coefficient distribu- 
tion corresponds to the expected random deviations 
for the Gaussian disturbances in the ACDM model, 
we may arrive at a conclusion that there exists an in- 
finitesimal number of separated non-Gaussian spots, 
if any at all. And allowing for the fact that the Cold 
Spot itself might be the result of a local modulation of 
spherical harmonics, that appeared as a consequence 
of the component separation ^33j, we may construe 
that this phenomenon, observed simultaneously in 
CMB and NVSS, is most likely not an anomaly but 
may be a systematism (foreground overestimation) 
while analysing the WMAP data and a random co- 
incidence with a statistically insignificant peculiarity 
of the NVSS distribution. 

5. CONCLUSION 

We have presented here a method of correlation map- 
ping on the sphere for a given angular scale. This 
method is software-implemented within the frame- 



work of the GLESP package. The proposed method 
allows analyzing on the sphere the properties of a 
random CMB signal having a single realisation based 
on its statistical characteristics only, namely, on its 
ergodicity, when taking a multitude of CMB realisa- 
tions in different regions of the sphere, we can make 
a conclusion about its realisation in a multitude of 
similar Universes, and, in so doing, evaluate its prob- 
able values. We have demonstrated the capabilities 
of our method for the WMAP5 data (the maps of 
ILC and foreground radiations: synchrotron, free-free 
and dust). For the correlation maps constructed, we 
discovered a shift in the pixel value distribution in 
favour of anti-correlations for dust and synchrotron 
radiation on the scales of 162', 324', and 540'. This 
calls for a possible overestimation of the contribution 
of these foregrounds in the method of ILC compo- 
nents separation. An excess in the negative coefficient 
values indicates that in the corresponding pixels, the 
correlated signals of CMB and dust have an oppo- 
site sign. This fact may be used while correcting the 
contribution of the dust component. We are planning 
further studies of this issue. 

The described approach was as well used for a 
study of correlation properties of the CMB and NVSS 
maps. In the correlation coefficient distribution there 
was found no additional signal a level higher than la 
over the one expected in the event of a random agree- 
ment of Gaussian disturbances in the ACDM cosmo- 
logical model. Such a distribution of correlations ar- 
gues for the fact that the Cold Spot, observed both 
in the CMB and NVSS, may be in fact the result of 
a random coincidence, and there is no need invoking 
new physics to explain it. 
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